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Using first-principles density functional calculations, we study the electronic and magnetic proper- 
ties of ferromagnetic insulating double-perovskite compound L^NiMnOe, which has been reported 
to exhibit interesting magnetic field sensitive dielectric anomaly as a function of temperature. Our 
study reveals existence of very soft infra-red active phonons that couple strongly with spins at the 
Ni and Mn sites through modification of the super-exchange interaction. We suggest that these 
modes are the origin for observed dielectric anomaly in L^NiMnOg. 

PACS numbers: 71.20.-b,75.30.Et,71.45.Gm 



Double perovskite L^NiMnOg (LNMO) is an inter- 
esting compound which is a ferromagnetic (FM) insu- 
lator with a Curie temperature close to room tempera- 
ture (Tc ~ 280 K). Recently large magnetic field induced 
changes in dielectric properties have been observed in 
this compoundi, which makes this compound a promis- 
ing material for potential device applications^. In spite 
of great technological importance, the theoretical effort 
for understanding this material is rather limited. To our 
knowledge, there exists only one report of linear aug- 
mented plane wave based basic electronic structure cal- 
culations of LNMO, by Matar et alA 

In this letter, we carried out first-principles density 
functional calculations to understand the ferromagnetic 
insulating behavior in this compound as well as to un- 
derstand the origin of dielectric anomaly that has been 
observed experimentally^. We used a combination of 
two types of methods, namely: (a) muffin-tin orbital 
based linear muffin-tin orbital 4 - and N-th order muffin- 
tin orbital (NMTO)^ methods, and (b) plane wave-based 
methods. In the latter, we used ultra-soft pscudopotcn- 
tials with an energy cutoff of 25 Ry (150 Ry) on the 
plane wave basis for wave functions (charge density) and 
a 6x6x6 mesh of k-points in sampling the Brillouin zone 
for a phase with unit cell containing one formula unit 
and equivalent for other phases. In particular, the struc- 
tural optimization and phonon calculations have been 
carried out using QUANTUM ESPRESSO^ and effective 
charges and dielectric response have been carried out us- 
ing ABINIT— . The analysis of hopping interactions by 
constructing effective orbitals, on the other hand, has 
been carried out within the framework of NMTO. In our 
LMTO and NMTO calculations, we have used four dif- 
ferent empty spheres to achieve the space filling. We 
used a spin polarized generalized gradient approximation 
(GGA)£ to the exchange correlation functional. 
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TABLE I: Energy-minimized structural parameters of LNMO. 
Lattice constants have been kept constant at the experimental 
values^. 

Rhomb ohedral 
a(A) b(A) c(A) x y z 

5.474 5.474 5.474 La 0.24980 0.24980 0.24980 

a [3 7 Ni 0.0 0.0 0.0 
60.671 60.671 60.671 Mn 0.5 0.5 0.5 

O 0.81403 0.67182 0.25889 
Monoclinic 

a(A) b(A) c(A) x y z 

5.467 5.510 7.751 La 0.00838 0.03781 0.24968 
a [3 7 Ni 0.0 0.5 0.0 
Mn 0.5 0.0 0.0 
90.000 90.119 90.000 01 0.22344 0.20903 0.04140 

02 0.29189 0.27756 0.45700 

03 0.42219 0.01454 0.24243 



LNMO, having the general structure of a double or- 
dered perovskite (A2BB'Oe), is distorted from the ideal 
double perovskite, and the amount of distortion changes 
as the temperature varies. The structure of La2NiMn06 
is rhombohedral (R3) at high temperature and trans- 
forms to monoclinic (P2i/n) at low temperature, with 
these two structures coexisting over a wide temperature 
ran g Q 9 i io j n v j ew f fh e f ac t t na t fh e positions of light 

atoms arc often not well characterized within the exper- 
imental technique, we have carried out structural op- 
timization of both rhombohedral (RH) and monoclinic 
(MC) phases where the internal degrees of freedom asso- 
ciated with La and O atoms have been optimized keep- 
ing the lattice parameters fixed at experimentally deter- 
mined values^. The relaxed structural parameters of the 
rhombohedral FM state (see Table I) agree well within 3 
% with experimental ones proving the reliability of our 
calculation scheme. The position of oxygen atoms, in 



particular, the x co-ordinate of 03 oxygen of the MC 
phase, however differ noticeably (a deviation of about 22 
%) from the experimental values (compared Table 3 in 
re^) . Our results may provide basis to further refinement 
of the experimental structure. Each NiOg octahedra in 
the rhombohedral phase is tilted with respect to MnO@ 
octahedra giving rise to Ni-O-Mn bond angle of 157°. 
The tilting is further increased by 2° in the monoclinic 
phase. We determined the electronic structure of geome- 
try optimized FM LNMO, for rhombohedral and mon- 
oclinic phases using the LMTO^ basis, as well as using 
the plane wave basis^. Both methods resulted in simi- 
lar features in the computed density of states and band 
structures, and insulating solutions for both RH and MC 
phases. The spin resolved partial density of states (DOS) 
of LNMO in the rhombohedral phase, calculated using 
ESPRESSO, is shown in upper panel of FigQ] Below - 
2 eV the predominant contribution is from 0-2p states. 
The octahedral surrounding of Mn and Ni atoms split the 
Mn and Ni d-manifolds into t2 ff and e g levels. In the up 
spin channel, the Ni-t2 S and Ni-e g levels are found in the 
energy range ~ -2 eV to Fermi energy and show a sig- 
nificant mixing with Mn-d states and 0-p states. In the 
down spin channel the Ni-t2 9 bands are located between 
0-p states and Fermi level, while Ni-e g states lye ~ 1.2 
eV above the Fermi level. This correspond to the nom- 
inal valence of Ni 2+ (d 8 : t^e 2 ,). In the up spin channel 
the Mn-t2 S bands are localized between Ni-t2 S and Ni-e s 
bands and arc filled, while the Mn-c 9 bands are sepa- 
rated by a gap of ~ 2.5 eV from the Mn-t2 9 bands and 
remains empty. In the down spin channel, both Mn-t2 9 
and Mn-e 9 bands are located above Fermi level in the 
energy range ~ 1.5 eV to 5 cV. This leads to conclu- 
sion that the oxidation state of Mn is nominally 4+ (d 3 : 
tlgGg), which agree with the Mn NMR and X-ray ab- 
sorption spectroscopy result s 12 ' 13 , though disagree with 
one of neutron diffraction studyJ^. Our spin-polarized 
LMTO calculations gave a moment of 3.0 fiB at the Mn 
site within a muffin-tin(MT) radius of 1.32 A, which agree 
with experimental value of 3.0 [Ie^- The magnetic mo- 
ment at the Ni site, for a MT radius of 1.52 Ais found to 
be 1.43 /xB, which is less than the experimentally mea- 
sured value of 1.9 fiB. The residual moment is found to 
reside at the oxygen sites giving rise to the total mag- 
netization of 5.0 ^B in agreement with the experimental 
valued. The spin resolved density of states in the mono- 
clinic structure is shown in the lower panel of Fig[TJ The 
basic nature of DOS is similar to that of rhombohedral 
phase of LNMO. The occupation of the Ni-d states and 
Mn-d states suggests again the nominal oxidation states 
of Ni and Mn ions to be 2+ and 4+, respectively The 
moments are found to be 2.91 fiB within a MT radius of 
1.38 Aat the Mn site, 1.35 /xB at the Ni site within a MT 
sphere of radius 1.52 A, and - 0.10 ^B (for 0.95 AMT 
radius) at the oxygen site, giving rise to ~ 5.0/f.u total 
magnetic moment, which is again in agreement with the 
experimental findings. 

LNMO being an insulator, the ferromagnetism in this 
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FIG. 1: GGA DOS of LNMO in geometry optimized rhom- 
bohedral and monoclinic phases. Zero of the energy is set at 
the GGA Fermi energy. 



compound is expected to be dominated by the localized 
super-exchange kind of interaction, resulting from the in- 
teraction of the half-filled d orbital of one metal ion with 
the vacant d orbital of another metal ion through an- 
ion p orbital. In the following, we show the feasibility 
of such a scenario by considering the calculated hopping 
interactions between effective Ni-d orbitals and Mn-d or- 
bitals. The construction of effective Ni- and Mn-d or- 
bitals have been carried out with NMTO- downfolding 
procedure, by integrating out O and La orbital degrees 
of freedom and keeping active only the Mn- and Ni-<i 
degrees of freedom. This procedure generates the ef- 
fective Ni- and Mn-d orbitals (see Fig. [2]) which takes 
into account the renormalization from the integrated-out 
oxygen and also La degrees of freedom. Considering an 
extended Kugel-Khomskii model^ that includes the hy- 
bridization between half-filled Ni e g orbitals and vacant 
Mn e g orbitals, the virtual hopping of parallely aligned 
spins is allowed and is favored over the virtual hopping 
of antiparallely aligned spins due to the energy gain via 
the Hund's coupling Jh- The net exchange can be ex- 
pressed in terms of the sum of the square of the Ni e g - 
Mn e g hopping interactions, J2 m m'(^e m e m ') 2 ' the ons he 
energy differences A ej6 , the onsite Coulomb U and the 
Hund's coupling Jh as: 



J (1) - -4 

J Ni-Mn ~ ^ 
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Considering the hoppings, and onsite energy differences 
for Ni-Mn neighbors, £ m m /(i e ™ em >) 2 turns out to be 
about 0.2 cV while A ej6 turns out to be about 1.9 eV. 
Considering a typical U value of 4 eV and Jh of 0.9 eV, 
it gives rise to a value of about -24.5 meV for the FM 
exchange interaction between Ni and Mn. 

However, the Ni-Mn exchange interaction has contri- 




FIG. 2: (Color on-line) Overlap between effective x 
orbitals, placed at neighboring NiOe and MnOe octahedra 
of LNMO calculated in the monoclinic phase, showing the 
super-exchange path mediated by the corner-shared oxygen. 
Plotted are the orbital shapes (constant-amplitude surfaces) 
with lobes of opposite signs colored as black(magenta) and 
white(cyan) respectively for Mn(Ni). 



bution also from the interaction between half-filled Ni e g 
orbitals and Mn t 2g orbitals, which should be antifcrro- 
magnctic (AFM) in nature, given by: 



y 
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The summation to, to' in the above runs over half-filled 
e g and t 2g orbitals at Ni and Mn sites, respectively. The 
computed sum of squares of the Ni e g - Mn t 2g hopping in- 
teractions in the basis of NMTO-Wannicr functions turns 
out to be about 0.02 cV, while the on-site energy differ- 
ence turns out to be about 0.25 eV. Putting these values 
in Eqn. (2) it gives rise to a value of about 19 meV for . 
The net exchange interaction, therefore, comes out to be 
ferromagnetic with a value of about 5 meV— , which is a 
reasonable estimate considering the rather approximate 
nature of the perturbative approach. The mean-field T c 
computed with this estimate of Ni-Mn exchange interac- 
tion comes out to be 350 K compared to experimental 
estimate of 280 K±£. 

In the following, we investigate the origin of mag- 
netocapacitance in LNMO manifested in the dielectric 
anomaly^ as a function of magnetic field. The dielec- 
tric constant of LNMO is known^ to increase with tem- 



perature and exhibit a jump at T,- 



T, 



depends 



sensitively on magnetic field, particularly for small fields 
(< 0.1 T). Since the electronic contribution to dielectric 
constant of insulators such as LNMO is typically much 
smaller— than the magnitude of jump in the dielectric 
constant, we expect the origin of this coupling between 
magnetic field and dielectric response to emerge from the 
couplings between spin and structure, i.e., phonons. 

To determine the coupling between spin and various 
phonons, we studied the response of optimized FM rhom- 
bohcdral structure to changes in magnetic ordering: e.g. 
changes in phonon frequencies with changes in magnetic 
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FIG. 3: (Color online) Phonon spectra of rhombohedral 
LNMO in FM (top panel) and FIM (bottom panel) states. 
The arrows show the shifting of dominant IR active phonon 
modes. The inset shows the displacement of atoms corre- 
sponding to the lowest frequency IR-active mode. The angle 
between the red dotted lines connecting the Ni(at the center 
of the oxygen-octahedra)-0 and O-Mn (rightmost corner of 
the cell) is affected by this phonon. 



ordering. We determined changes in phonon frequen- 
cies upon changing the neighboring spin to antiparallel, 
the so called Ferrimagnetic (FIM) alignment. Note that 
any other different spin arrangement other than FM spin 
arrangement would have been equally qualified for this 
purpose. The charge states of Ni and Mn in the FIM 
state are found to remain same as that in the FM state. 
Hellman-Feynman forces acting on atoms in the FIM 
phase, give us the lowest order coupling between spins 
and phonons which is linear in atomic displacements. We 
find that these forces are equal and opposite for pairs of 
atoms, hence this coupling is zero for any IR-active modes 
and should have no direct implications to the observed 
dielectric anomaly. Next, we determined the T— point 
phonons for the rhombohedral structure with FM and 
FIM ordering. Shifts in the phonon frequencies give the 
coupling between spins and atomic displacements at the 
second order. Since the rhombohedral structure is un- 
stable at T=0 K, we find two marginally unstable modes 
(31i and 19i cm -1 ) which are IR-inactive and couple 
strongly with spins: their frequencies change to 64i and 
63i cm -1 in the FIM phase. We find that frequencies 
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of the lowest energy IR-active phonons soften from 91.3, 
149, 228 and 255 cm" 1 in the FM phase to 65.5, 120, 184 
and 199 cm -1 respectively (indicated by arrows in Figl3j), 
exhibiting a strong coupling with spin. This may be com- 
pared with the recently studied case of CdC^S/p^ where 
the significant polar mode was found only at frequency 
of 300 cm -1 . This results in change in static dielectric 
constant from 119 in FM state to 221 in the FIM state, 
which is dominated by the softest IR-active mode (see 
inset of Fig[3]) with a contribution of 77 and 185, respec- 
tively. Atomic displacements in this softest mode, are 
such that they would make the Ni-O-Mn angles closer 
to 180°, in an average sense, leading to enhancement of 
the supcrcxchangc interaction. The decrease in phonon 
frequency for the spin-paired Ni-Mn in FIM phase com- 
pared to spin-antipaired Ni-Mn in FM phase can be ex- 
plained by analyzing the spin Hamiltonian JS j .Sj and 
noting that the magnetic super-exchange coupling J de- 
pends on Z Ni-O-Mn, 9, as cos 2 (9). Expanding cos 2 (9), 
about the equilibrium value of Z Ni-O-Mn, 9q « 157°, 
and assuming 9 = 9q + u, u being the displacement, 
the spin-phonon coupling turns out to be positive for the 
term linear in u and negative for the term quadratic in u. 
The latter effectively gives a positive change in phonon 
frequency due to the additional negative sign, associated 
with FM nature of the magnetic interaction (cf. Eqn.l.). 

To conclude, we carried out first principles DFT calcu- 
lations to examine the electronic and magnetic structure 



of La2NiMn06, in particular, to understand the origin 
of the dielectric anomaly reported recently. We could 
correctly reproduce the ferromagnetic insulating behav- 
ior of the compound, the ferromagnetism being governed 
by the super-exchange interaction. Our study further 
showed presence of soft IR-active phonon modes exhibit- 
ing strong coupling with spin, which explains the observed 
dielectric anomaly. The fact that the jump in the di- 
electric constant at H=0 T, occurs at a temperature 
Tj ump below T c happens because close to T c the magnetic 
moment is not fully developed due to thermal fluctua- 
tion while at lower temperature the moment gets fully 
developed and makes the effect of coupling to phonon 
degrees of freedom appreciable enough to observe the 
jump in dielectric constant. This is corroborated by the 
fact that Tj ump becomes closer to T c upon application 
of magnetic field, which helps overcoming the thermal 
fluctuation and enhance the magnetization. Note that 
superexchange driven B-site magnetism based dielectric 
anomaly discussed here is complementary to the mecha- 
nism of the low-temperature dielectric anomaly discussed 
for EuTiO^, which is A-site based weak magnetism 
driven. 
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